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ABSTRACT 



Aims. We looked for signatures of Quasi-Biennial Periodicity (QBP) over different phases of solar cycle by means of acoustic modes 
of oscillation. Low-degree p-mode frequencies are shown to be sensitive to changes in magnetic activity due to the global dynamo. 
Recently have been reported evidences in favor of two-year variations in p-mode frequencies. 

Methods. Long high-quality helioseismic data are provided by BiSON (Birmingham Solar Oscillation Network), GONG (Global 
Oscillation Network Group), GOLF (Global Oscillation at Low Frequency) and VIRGO (Variability of Solar IRradiance and Gravity 
Oscillation) instruments. We determined the solar cycle changes in p-mode frequencies for spherical degree i=0, 1, 2 with their 
azimuthal components in the frequency range 2.5 mHz < v < 3.5 mHz. 

Results. We found signatures of QBP at all levels of solar activity in the modes more sensitive to higher latitudes. The signal strength 
increases with latitude and the equatorial component seems also to be modulated by the 1 1-year envelope. 

Conclusions. The persistent nature of the seismic QBP is not observed in the surface activity indices, where mid-term variations are 
found only time to time and mainly over periods of high activity. This feature together with the latitudinal dependence provides more 
evidences in favor of a mechanism almost independent and different from the one that brings up to the surface the active regions. 
Therefore, these findings can be used to provide more constraints on dynamo models that consider a further cyclic component on top 
of the 1 1 -year cycle. 

Key words. Methods: data analysis - Sun: helioseismology - Sun: activity 



1. Introduction 

Evidences in favor of shorter periodicities than the 1 1 -year pe- 
riod, ranging from days to years in solar activity proxies, are 
very well known in literature. For short-terms, the Sun often ex- 
hibits 27-day and 13.5-day periodicities, while the regime be- 
tween days and 1 1 years are referred as mid-term periodicities 
(iBaiL 120031) . Among these, particularly p rominent are the ones 
between 1.5-4 years dVecchio et aUl2009l) . We will refer to them 
as Quasi-Biennial Periodicity (QBP). This phenomenon ha s 
been observed in radi o flux on 10 . 7cm ( Belmont et all Il996l) . 



total solar irradiance ( Penza et all 



2006), flares and sunspots 



(lAkiokal [19871: iMursula et all 120031: iValdes-Galicia & Velasco . 
' 20081) and green coronal emission line (IVecchio & CarboneL 
20081: IVecchio et all 120091) . It also seems to appear only time 
to time tending to be more important over periods coinciding 
with solar max ima and are not recognized as typical feature 
of every cycle dKriyova & SolankiL 120021: iKnaack et all 12004 
ICadavid et aILl2005HLi et al.Ll2006l) . This suggests that the QBP 
might be closely linked to the strength of the global dynamo. 
Coronal dat a have shown a latitudinal dependence of the sig- 
nal strength dVecchio & Carbonell2008l) . It is particularly strong 



in the polar region and weaker around the equator and more 
over the periodicity from northern and equatorial regions (s;2.8 
years) is different from the one in southern regions («1.5 years). 
More recently the same authors, investigating a longer dataset 
recorded from 1939 to 2005, concluded that the non-constant 
period length s is the manifestation of a unique quasi-biennial 
activity cycle ( IVecchio et all 120091) . 

Since the discovery of the correlation between the observed 
global acoustic mode (p- mode) parameters a nd the solar mag- 



netic activity cycle (e.g. Woodard & Novesl 11985b iPalle et all 
119891 : lElsworth et all 1 19931) . it became evident that the acoustic 
modes can also be used as a diagnostic tool of the solar cycle. In 
particular with high-quality and continuous helioseismic data, 
deeper analysis revealed that mode amplitudes are suppressed 
with increasing magnetic activity, while mode frequencies are 



maximum at low and high-angular degree (Howe etall 1999 


Kommetal., 20001 


Jimenez-Reyes et al., 2001; Chaplin et al. 


20011: ISalabert et al 




2004: Salabert & Jimenez-RevesL 2006 


Simoniello etall 120 





). While for the mode amplitudes a possi- 



ble explanation in terms of mode conversion has been provided 
dSimoniello et all 1201 Ol) . we still lack the basic understanding 
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of the mode frequency shift. F-modes can help to understand 
this variability and the confrontation of obse rvations to pre- 
dictio ns seem to show some structural effects dLefebvre et al.L 
120091) . Unfortunately, a proper introduction of sub-surface mag- 
netic field is still lacking. So this does not allow to properly 
estimate the structural consequences of this important factor. 
The form and the degree dependence of the shifts, anyway fa- 
vor near-subsurface phenomenon, but however they cannot be 
purely explained by structural chang es, the magnetic field is 
some h ow involved in the mechanism dGoode & DziembowskiL 
l2002t iDziembowski & Goodel |2004|) . Recently the argument 
that the solar cycle variation in the global mod e frequency are 
due t o the global averaging of active regions (iHindman et al.L 
1200 lb has been only partially supported by the observations 
( Tripathv et al.L 120101) . It has been argued that the weak com- 
ponent of the magnetic field (turbulent or horizontal field) must 
be taken into account in order to fully explain the observed mode 
frequency shift. 

Recent observations in p-mode frequencies have also shown 
the presence of a 2-year modulation most prominent over peri- 
ods coinciding with solar maxima. The properties of the seis- 
mic QBP have been interpret ed as the visible manifestation of 
a second dynamo mechanism dBenevolenskaial I1995L Il998alfbh . 
induced by the strong latit udinal shear locate d below the sur- 
face at su 0.95 solar radius (ISchou et al.L 1 1998b . The second dy- 
namo itself cannot explain the observed modulation in the signal 
strength and therefore it has been supposed that when the 11- 
year cycle is in its strong phase, buoyant magnetic flux are sent 
upward by the main dynamo into shallower layers, thus allow- 
ing the s eismic QBP to be ma inly detected over periods of solar 
maxima dFletcheretalll2010b . However it is still a matter of de- 
bate if this second dynamo mechanism is indeed behind the QBP 
observed in activity proxies and in solar seismic data. In fact, re- 
cent observations of rapidly rotator stars have shown signatures 
of QBP, whose strength increases with increasing magnetic ac- 
tivity and it gets stronger at polar regions. It has been supposed 
that the unst able magnetic Rossby waves in t he tachocline might 
be behind it dZaqarashvili et al.L 120101 [201 ll) . If this is also true 
for solar type star is still an undergoing investigation. 

The p-mode spatial configuration can be described by spher- 
ical harmonics with each mode being characterized by its spher- 
ical harmonic degree { and azimuthal order m. Modes with I < 3 
are named low-degree modes and conversely to surface activity 
indices, they also sound the interior of the Sun being shaped in 
deeper layers. Therefore they might feel magnetic changes that 
have yet to manifest at the surface or it might also be possible 
that the disturbances will never reach the Sun's surface. By av- 
eraging the individual p-mode frequency shifts over all low {€)- 
values, observational findings pointed to seismic signatures of 
QBP over the long extended minimum charact erizing the end of 
solar cycle 23 and beginning of solar cycle 24 dBroomhall et al.l 
l2009aL iFletcher et all L20Toh . Further evidences in favor of this 
feature have been reported in low-deg ree modes by investigat - 
ing the spherical degree dependence dBroomhall et all 1201 ll) . 
but some authors have cast doubts, because they didn't spot 
any s ignificant trend during periods of low activity dJain et al.L 
[201 lb . We decided, therefore, to investigate the properties of the 
seismic QBP by using several helioseismic instruments in veloc- 
ity and intensity data covering a period that goes from the end 
of solar cycle 22 up to the beginning of solar cycle 24. Each 
observational program can provide more sta ble and/or accu- 
rate data depending on instr umental features dBroomhall et al.L 
l2009btlSaiabert et al.L [20 1 lb . Therefore by comparing different 
observations we have been able to check whether or not the seis- 



mic QBP occurs only time to time and during periods of high 
activity. Furthermore we carried out this analysis for each of 
the sectoral (m-{) and zonal (m=0) components of low-degree 
modes (without averaging over the azimuthal components) and 
this is the novelty of the work within this topic. We also investi- 
gated if and to which extent the mid-term periodicity is linked to 
the appearance of active regions induced by the main dynamo. 
To this aim we carried out a correlation analysis between each 
azimuthal (m)-components of low-degree mode frequencies and 
activity indices linked to strong and/or weak magnetic fields over 
different phases of solar cycle 23. It is important to investigate 
the phasewise variation in correlation between low-{ frequen- 
cies and activity indices as well, because earlier study of the m- 
componen t s were confined to intermediate degree modes only 
dJain et al.L 12009b . Furthermore, acoustic waves, depending on 
their spherical degree and azimuthal component, are more sensi- 
tive to low-mid latitudes or higher latitudes. Therefore the anal- 
ysis of the azimuthal components allowed us also to have some 
hints on the role played by strong and/or weak fields with lati- 
tudes as sources of the QBP. 

The seismic investigation of the mid-term periodicity might 
help us to gain a deeper insight on the global properties of the 
Sun as for example on the solar dynamo theory as well as on stel- 
lar cycles: mid-term periodicities have already been observed in 
fast rotator stars and therefore the better understanding of this 
short periodic behavior in the Sun might help illuminate the in- 
terpretation of this signal of solar-like stars. 

2. Helioseismic data analysis 

2. 1 . Mode visibility 

Low-degree modes are detected either in Sun-as-a-star observa- 
tions or through highly spatially resolved images. In the first sce- 
nario the modes are observed in averages made over the visible 
disk of Doppler velocity or intensity variations. In this case only 
the modes that satisfy the condition 6 + m = even are visible, 
because of the near-perpendicular inclination of the Sun with re- 
spect to the Earth. In observations with high spatial resolution, 
the visible disk is sampled in many pixels and the collected im- 
ages are decomposed into their constituent spherical harmonics, 
giving access to all 21+ 1 component for each degree. 

2.2. Sun-as-a-star observations from ground based network 

Birmingham Solar Oscillation Network (BiSON) has been op- 
erating since 1976 and makes unresolved solar disc observa- 
tions. It consists of 6 resonant scattering spectrometers, that per- 
form Doppler velocity measurements in integrated sun light of 
the K Fraunhofer line at 7699 A dChaplin et al.L 1 1995b . The data 
are dominated by the Doppler variations from the low-degree, 
(, modes. T he nominal height is at ^26 km above the photo- 
sphere (e.g. iJimenez-Reves et ail 12007b . We analyzed the data 
provided by BiSON instruments from 1 1th April 1996 to 4th of 
April 2009. This long dataset was split into contiguous 365-day 
series with four overlaps of 91.25 days. Estimates of the p-mode 
frequencies for ^=0,1,2 were extracted from each subset by fit- 
ting a modified Lorentzian mode l to the data using a standard 
likelihood maximization method (Fle tcher et all 12009b . 

2.3. Sun-as-a-star observations from space 

The Global Oscillation at Low Frequency (GOLF) instrument, 
on board the SOlar and Heliospheric Observatory (SOHO) satel- 
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lite, is devoted to the search of low-degree modes (Gabriel et 
al. 1995). It works by measuring the Doppler shifts of the Na 
line at 5889A (Dl) and 5896A (D2). Due to the malfunctioning 
of the polarization system, GOLF is working in a single-wing 
configuratio n and it has been p roved that it is almost a pure ve- 
locity signal dPalle et al.Lll999h . The one-wing configuration has 
been changing during the 14 years of observations as follows: 1) 
1996-1998 in the blue-wing configuration; 2) 1998-2002 in the 
red- wing configuration; 3) up today in the b lue-wing configura- 
tion (lUlrich et all l2000L iGarcfa et alll2005h . 

Due to the formation heights of the spectral lines, in the 
blue-wing configuration, GOLF observes at ^330 km, while in 
the re d-wing configuration at ^480 km dJimenez-Reves et"ail 
l2007h . Although the changing o bservational height affects the 
p-mode amplitude dete rmination ( ISimoniello & Salabertl 120091 : 
ISimoniello et all [201 Oh . it doesn't affect the p-mode frequency. 
We analyzed the data provided by GOLF instruments from 1 1th 
April 1996 to 4th of April 2009 and the estimates of p-mode 
frequency shifts have been extracted by applying the same fit- 
ting method as in BiSON. This differs from the one applied to 
VIRGO and GONG £-0 integrated. We, therefore, verified that 
both fittings methods returned consistent results. 

The Solar PhotoMeters (SPM) of the Variability of so- 
lar IRradiance and Gravi ty Oscillations (VIRGO) package 
dFrohlich et all 1 1 995L fl997l) aboard SOHO satellite, consists of 
three independent photometers, centered around 402, 500 and 
862 nm (the blue, green, and red channel respectively). They 
measure the spatially integrated solar intensity over a 5 nm 
bandpass at a one minute cadence. The data obtained by the 
VIRGO/SPM instrument over the mission have been of uni- 
formly high quality, whether before or after the loss of contact 
with the spacecraft for several months in 1998. Each one of the 
three different signals observes the Sun- as-a-star oscillations at 
different heights in the solar atmosphere (I Jimenez et alll2005h . 

The data provided by VIRGO instruments cover almost 
14 years of data, starting on the 11th April 1996 up to 8th 
April 2010. This dataset was split into contiguous 365-day 
with four overlaps of 91.25 days. The power spectrum of each 
subseries was fitted to estimate the mode p-mode parameters 
for £=0,1,2 using a s tandard likelihood maximization function 
dSalabert et all 120071) . Each mode component is parameterized 
using a asymmetric Lorentzian profile. 

2.4. High-resolved observations from ground based network 

The Global Oscillation Network Group (GONG) consist of 
six instruments deployed worldwide to provide nearly contin- 
uous and stable velocity images of the Sun. GONG instru- 
ment is base d on a Michelson i nterfe rometer called a Fourier 
Tachometer dBeckers & Browrl QJ)78). It works by using the 
Ni line at 676. 8nm and the observational height is approx- 
imately at 213 km. The network started to be fully operat- 
ing since May 1995. Time series of 36 days (the so-called 
GONG month) are produced and the individual frequencies 
of the azimuthal components are determined up to I = 150 
over subseries of 108 days by concatenating three GONG 
months. These frequencies are afterwards made publicly avail- 
able ( http : //gong2 .nso. edu/archive/) . 

The i=0 GONG time series between May 7th 1995 and May 
3rd 2008 was analyzed in a similar manner as for VIRGO data, 
returning estimates of the mode frequencies of ^=0,1,2. In order 
to analyze the azimuthal components individually, we made use 
of the frequency tables of the (=1,2 modes available from the 
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Fig. 1. Solar cycle changes in p-mode frequency shifts (pHz) for 
£=0, 1, 2 from VIRGO (green), BiSON (black), GONG (pur- 
ple) and GOLF (blue and red to identify the corresponding wing 
configuration) observations over solar cycle 23 and beginning of 
solar cycle 24. 



GONG webpage between the 29th of June 1995 and the 24th 
of December 2009. Very few outliers were removed from the 
following analysis using the quality flags provided by the GONG 
team. 

2.5. Frequency shift determination 

The temporal variations of the p-mode frequencies were defined 
as the difference between the frequencies of the corresponding 
modes observed at different dates and reference values taken as 
the average over the years 1996-1997. The weighted averages of 
these frequency shifts were then calculated between 2500 and 
3500 fiHz. Furthermore for each spherical degree the shifts from 
the VIRGO green, blue and red channels have been averaged. 
The shifts for each spherical degree with the same | m | obtained 
by the azimuthal analysis of GONG high-resolved data have also 
been averaged. 

3. QBP signatures in p-mode frequencies over solar 
cycle 23 

3. 1 . Spherical degree dependence of the QBP over the 
11 -year cycle 

We looked for the existence of QBP signatures in p-mode fre- 
quency shifts over the 11 -year cycle. Fig. [1] shows the spherical 
degree dependence of the shift as seen from VIRGO, BiSON, 
GONG and GOLF datasets. We introduced an artificial offset 
between the curves to better compare the seismic features. We 
can spot several seismic signatures that might be linked to the 
QBP signal: 

- for £=0 modes over the minimum of solar cycle 22 and 
beginning of solar cycle 23 (1995-1997), GONG observations 
show a seismic signature of a QBP. Following the ascending 
phase of solar cycle 23 in BiSON and GONG, between 1997- 
1999, we find a further modulation of about 2 years. This signa- 
ture doesn't appear in GOLF and VIRGO data, probably due to 
the loss of SOHO in June 1998. During the maximum of solar 
activity the modes show up with a clear double peak structure 
in all datasets. Over the descending phase BiSON, GOLF and 
VIRGO agree quite well, showing the same downward trend not 



3 



R. Simoniello et al,; QBP from helioseismic observations 



characterized by any modulation. GONG data, instead, shows 
an oscillatory signal around December 2003. Towards the end of 
solar cycle 23 and beginning of solar cycle 24 (2006-2009), we 
observe a QBP in all dataset although in GOLF and VIRGO data 
are particularly prominent; 

- in i-\ modes there are no traces of QBP over the ascending 
phase in BiSON, GOLF and VIRGO data. Instead GONG shows 
a modulation around July 1998. Over the maximum phase we 
found a clear double peak structure in all dataset, but in GONG 
is more pronounced. In GONG, BiSON and GOLF, over the 
descending phase, we might spot two further modulations be- 
fore December 2003 and before September 2006. that are in- 
stead absent in VIRGO data. Over the long extended minimum 
of solar cycle 23 only VIRGO shows a further periodicity after 
September 2006; 

- in 1-2 modes we found over the minimum of solar cycle 22 
and beginning of solar cycle 23 (1995-1997) a seismic signature 
of a QBP from GONG observations. We found the same sig- 
nature in GONG data in t-Q mode. Over the ascending phase, 
a strong peak of almost 2-year length around 1998 in VIRGO 
and GOLF data is present. This signal might be an artefact due 
to the loss of the satellite between June and September 1998. 
Over the maximum phase we found a double peak structure in 
BiSON, GONG and GOLF observations. The descending phase 
is characterized by a very well modulated signal of about 2 years, 
although in BiSON, GOLF and VIRGO is more pronounced. 
During the minimum of solar cycle 23, we spot in BiSON and 
GOLF data a clear signature of a further periodicity starting in 
September 2006 and ending before June 2009. 

In summary we found signatures of biennial modulation in 
all observational programme, but there are differences in the data 
sets. The signal strength of mid-term variations is weaker com- 
pared to the 1 1 -year signal. In a recent paper the mode visibility 
in GOLF data appears different in the two wing- configurations , 
being slightly better during red-wing period (ISalabert et all 
1201 lb . This might explain why the faint signal of the QBP is 
less prominent in GOLF data. For intensity observations further 
aspects need to be taken into account, like the selection of the 
spectral region in the observations that affect the mode v i sibilit y 
of the modes in photometric measurements dNutto et all l2008h . 
Furthermore the continuum intensity data samples a different 
(i.e. deeper) region where the temperature variations induced by 
the acoustic oscillations may be quite small. 

By analyzing the spherical degree dependence of p-mode 
frequency shift, we detected likely signatures of the 2-year peri- 
odicity at all levels of solar activity, instead surface activity in- 
dices mainly show it over periods coinciding with solar maxima. 
We also noticed that away from the solar maximum, the QBP 
shows up in each low degree mode over different periods of so- 
lar activity phase. This results from the different spatial config- 
uration of the modes over the solar disk. For example, we found 
seismic signatures over the low activity phase in the modes (-Q, 
2 but not in the dipole mode I- 1 . In integrated sunlight observa- 
tions due to the mode visibility, the observed l-\ frequency is a 
weighted measurement of the visible components ((=1, \m\=l), 
while in the case of the 1-2 mode, the fitted frequency corre- 
sponds to the weighted measurement of the zonal {(-2, m-Q) 
and s ectoral ( (=2, \m\=2) components dJimenez-Reves et"al 



120041) . The spatial structure of the oscillation means that the 
sectoral \m\=t are more concentrated around equatorial latitudes 
than the zonal components (m=0), and is pa rticularly true with 
increasing ( (IChristensen-Dalsgaardl 120021) . As consequence 
the contribution to the QBP signal over the minimum activit y 
might have its origin at higher latitudes (ISalabert et all l2009h . 
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Fig. 2. Solar cycle changes in p-mode frequency shifts (yiHz) for 
the m-components of (-1, 2 modes as seen from GONG obser- 
vations. 



Furthermore the presence of the QBP signal over periods of low 
activity is not peculiar of so lar cycle 23, as guessed by some 
authors dFletcheret all 1201 Oh . 



3.2. Azimuthal dependence of the QBP over the 1 1 -year 
cycle 

We used GONG data to look for QBP signal in the azimuthal 
components of the dipole (£=1) and quadrupole {1-2) mode. 
Fig.|2]shows the solar cycle changes in /5-mode frequency shifts. 
The sectoral components exhibit a stronger 11 -year envelope 
compared to the zonal ones. This finding clearly confirms the 
different sensitivity of the modes to magnetic changes with lati- 
tude. The £=\, m-Q seems to show a mid-term periodicity over 
the end of solar cycle 23 with an upward trend starting soon after 
September 2006. If we look, instead, at t-\ m=l component, we 
spot the double peak structure over the maximum of solar activ- 
ity as also observed in integrated sunlight observations. A further 
mid-term periodicity seems to occur across September 2006, and 
the same structure appears in GONG (integrated) and BiSON 
data. Looking at both m-components of the quadrupole mode, 
the solar maxima is also characterized by the typical double peak 
structure. Then over the descending phase, the m—2 component 
shows one smoother signature of QBP across September 2006 
as in (=\, m=l. The zonal component, instead, seems to show 
a further periodicity starting soon after September 2006 as we 
already found in the mode (—1, m-Q. Therefore the azimuthal 
analysis brings more evidences that over the minimum activity 
phase of solar cycle 23 the QBP might have occurred at high 
latitudes. 



4. Latitudinal dependence of the QBP signal 
strength 

4.1. Modulation in the QBP signal strength 

To better reveal the seismic signatures of the QBP signal, we 
took out the 11 -year d ependence by using a boxcar of 2.5 
years, as already done in dFletcher et alll201QI) . Fig.[3]shows the 
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Fig. 3. The QBP signal in each low-degree p-mode frequency 
shifts after the dominant 11 -yr signal has been removed. The 
color legend is the same as Fig.Q] 



changes in p-mode frequency shifts due to the presence of the 
QBP signal for ^=0, 1, 2 over solar cycle 23. GONG, BiSON, 
VIRGO and GOLF observations point to a modulation of the 
signal strength in phase with the solar cycle in agreement wit h 
previous findings dFletcher et al l l2010tlBroomhall et al.Ll20TTb . 
It also seems that the modulation affects more the {=1,2 modes 
compared to i=0. This is a further evidence of the higher sen- 
sitivity to low-mid latitudes of the sectoral modes. Therefore 
we might also argue that the observed QBP signal in the sec- 
toral modes might come mainly from low-mid latitudes. If it is 
indeed the case, the signal at equatorial regions is of compa- 
rable strength (or even larger) to the QBP signal produced at 
high latitudes. This conclusion won't be correct at this stage. 
The magnetic field, depe nding on its strength and inclination 
(ISchunker & CallvLl2006h . affects the size of the acoustic cavity, 
where the mode propagates. But as the acoustic cavity decreases 
when the spherical degree increases, therefore, to get informa- 
tion on the strength of the signal coming at different latitudes, 
we have to compare sectoral and azimuthal components of the 
same spherical degree. Fig. [4] shows the QBP signal over so- 
lar cycle 23 for the m-components of the dipole and quadrupole 
modes. The residuals for 6=1, m=Q might indicate the presence 
of a QBP across July 1998 and soon after September 2006. If we 
look, instead, at 6=1 m=l component, we spot a nice-increasing 
trend in the signal strength of the QBP over the ascending and 
maximum phase of solar cycle 23. A further periodicity seems to 
occur across September 2006, and the same structure appears in 
GONG (integrated) and BiSON data. Looking at the zonal com- 
ponent of 6=2, we find over the solar maxima the typical dou- 
ble peak structure and a further mid-term periodicity soon after 
September 2006 as in €=1, m=0. The residuals for the sectoral 
component shows a nice-increasing trend in the signal strength 
of the QBP over the ascending phase of solar cycle 23. A further 
2-year modulation might have been occurred across September 
2006 as found also in 6=1, m=l. The azimuthal analysis con- 
firms the modulation in the signal strength mainly occurring in 
the sectoral components of the modes and also that over the long 
extended minimum of solar cycle 23 the sources of the mid-term 
periodicity might be located at high latitudes. 
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Fig. 4. The QBP signal in the zonal and sectoral component of 
the dipole and quadrupole mode frequency shifts after the dom- 
inant 1 1-yr signal has been removed from GONG high-resolved 
observations. 



4.2. QBP signal strength over high and low activity phase 

The zonal and sectoral components of the quadrupole mode 
show clear signatures of the QBP signal over periods coincid- 
ing with solar maxima. By comparing the peak-amplitude of the 
ra-components over this activity phase, the signal seems to be 
of comparable strength at different latitudes. But as the strong 
magnetic fields are confined at low-mid latitudes, the compari- 
son should be carried out over periods of low activity, when the 
magnetism of the Sun at low as well as high latitudes is dom- 
inated by weak fields. To this aim we might attempt to com- 
pare the strength of the possible signatures of the QBP signal we 
found over the years 2006-2009. For the quadrupole mode the 
strength is larger for the m=Q compared to m=2 component and 
the same feature occurs in the dipole mode. As consequence this 
might be a further hint that the QBP signal at high latitudes is 
larger than the one at low-mid latitudes. 

4.3. Investigating the significance of the QBP signatures over 
solar cycle 

Evidences in favor of QBP have been found in all the 
modes. In order to investigate their significance, we ap- 
plied the wavelet analysis developed by Torrence and Compo 
dTorrence & Compel 1 19981) to the spherical and azimuthal com- 
ponents of the modes. We show the results only for GONG ob- 
servations, as the findings for the spherical degree analysis from 
the different datasets are in good agreement. Fig. [5] shows the 
wavelet analysis carried out for each spherical degree. Within 
the cone of influence the local wavelet power spectrum identi- 
fies the QBP signal over the years 1998-2004 of solar cycle 23 
in all low-degree modes. The significant periodicity at 90% of 
confidence level occurs at T=l .7 for ^=0,2, while for 6=1 we 
have two peaks above the 90% confidence level corresponding 
to T=1.7 and T=2.1. Fig. [6] -[7] shows the wavelet analysis for 
the azimuthal components of the modes. In this case we dropped 
down the significant level at 80% to show that there are signa- 
tures of QBP signal over the minimum activity phase in the zonal 
components of the modes. The significant periodicity at 80% of 
confidence level occurs at T=l .7 for €=l,m=Q, T=2.1 for €=1,2 
and m=l,2 and T=1.6 for €=2,m=Q. It is interesting to note that 
the zonal modes show a slightly different periodicity from the 
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1995 1999 2003 2008 0.000 0.001 0.002 1995 1999 2003 2008 0.000 0.001 0.002 1995 1999 2003 2008 0.000 0.005 0.010 
Time Power Time Power Time Power 

(a) (b) (o) (b) (a) (b) 

Fig. 5. (a) The local wavelet power spectrum for (-0, 1, 2. White represent areas of little power, red those with the largest power. 
Solid rings identify areas at power of 0.001, 0.002,0.005 over time and only those that achieved 90% confidence level are labelled. 
(b)Global wavelet power spectrum. The dotted line are the 90% of confidence level. 




1995 1999 2004 2009 0.000 0.005 0.010 1995 1999 2004 2009 0.00 0.01 0.02 
Time Power Time Power 

(o) (b) (o) (b) 

Fig. 6. The same as in Fig. 5 for the t- \ m=Q,l. Solid rings identify areas at power of 0.0050,0.0075,0.01 over time and only those 
that achieved 80% confidence level are labelled. 




one detected in the sectoral components. This might be the result 
of the greater sensitivity to higher latitudes of the zonal modes. 
To further investigate this point, it will be important to study 
the properties of the QBP signal in the azimuthal components 
of intermediate-high degree modes. There are evidences, in fact, 
showing that the periodicity from northern (50° < 9 < 90°) and 
equatorial (-45° < 6 < 45°) region s differs from the one in th e 
southern regions (-50° < 6 < -90 0s ) dVecchio & CarboneLl2008l) . 

In summary the wavelet analysis confirms that a further pe- 
riodicity of k 2-years is indeed present in the data at confidence 
level of 90%. Furthermore the findings seem also to point to a 
likely persistent nature of the QBP signal, since we found evi- 
dences that this mid-term periodicity is also present over periods 
of low activity phase. 



5. p-mode frequency shift correlation coefficient 
with activity proxies over the 11 -year cycle 

5.1. Spherical degree dependence of the correlation 
coefficient 

The correlation analysis between p-mode frequency shift and 
surface activity measures might improve our understanding on 
the role played by strong and weak fields in the QBP signal. 
We, therefore,decided to investigate the degree of correlation be- 
tween /9-mode frequency shift and two activity proxies: MWSI 
(Mount Wilson Sunspot Index) and F10.7 index. The MWSI val- 
ues are determined as a summation over pixels where the ab- 
solute value of the magnetic field strength is greater than 100 
gauss. We have chosen MWSI because it is linked to the strong 
component of the toroidal fields generated by the global dynamo. 
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Fig. 8. The degree of correlation between p-mode frequency 
shifts of low degree modes (^=0,1,2 from left to right) and 
MWSI (top panel) and RF (bottom panel) over different phases 
of solar activity. The colors bar identify the observational pro- 
gramme, whose legend is the same as Fig.Q] 



These fields are mainly located at low-mid latitudes. F10.7 (RF) 
represents, instead, an integrated activity proxy and tells us the 
irradiance at 2.8 GHz at chromospheric heights coming mainly 
from the Quiet Sun background emission and also Sunspots, 
Plages. This proxy therefore is a good measure of the weak (lo- 
cated at high latitudes) and strong magnetic fields (mainly at 
low-mid latitudes). 

We investigated the degree of correlation of p-mode fre- 
quency shift with MWSI and RF by splitting the whole cycle 
in different phases as follows: 

- AP= Ascending Phase from 22 September 1996 until 26 
June 1999 

- MP= Maximum Phase 1999 June 27 to 2003 January 12 

- DP= Descending Phase January 13 2003 26 until July 2007 

- mP= minimum Phase up to 2009. 

We used the Spearman's rank correlation in preference to 
the Pearson's correlation because the Pearson's rank correlation 
assesses the linear correlation between two sets of data. The 
Spearman's rank correlation is less specific. While some of the 
data sets we compare may be linearly related this is not defi- 
nitely true in every case. Fig.[8]shows the correlation coefficient 
for €—0, 1, 2 mode frequency shifts with MWSI (upper panel) 
and RF (bottom panel) over different phases of solar cycle for 
BiSON, GONG, GOLF and VIRGO. We find that the degree 
of correlation of p-mode frequency shift with MWSI decreases 
over the maximum phase independently from I value. If we look 
at RF we find higher correlation compared to MWSI all over dif- 
ferent phases of the sol ar cycle. These fi ndings are in agreement 
with previous results dJain et all 120091) and they confirm that 
behind the mode frequency shift the active regions are not the 
only players. The other important feature is the anticorrelation 
over the minimum of solar cycle 23 for £=0, 2. In these modes 
there are traces of the QBP signal. This might be the result of the 
different behavior of the surface activity measure over the long 
extended minimum. While MWSI went down to unprecedented 
level until 2009, Radio Flux shows already an upward trend in 
February 2008. This might explain the higher degree of anticor- 
relation with RF compared to MWSI. 
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Fig. 9. The degree of correlation between p-mode frequency 
shifts of the azimuthal components of the modes and MWSI and 
RF. The colors bar identify MWSI (pink) and RF (purple). 



5.2. Azimuthal dependence of the correlation coefficient 

It is possible to investigate the latitudinal dependence of the 
correlation coefficient by using the single m-components of the 
modes provided by GONG data. 

Fig.|9]compares the degree of correlation between the zonal 
and sectoral components of dipole and quadrupole modes with 
MWSI and RF. The sectoral components seems to correlate bet- 
ter with both activity proxies over different phases of the so- 
lar cycle compared to the zonal modes. Previous works have 
shown that the central frequency of each multiplet is best cor- 
related with the global activity measure, while the individual 
m-component frequenci es are more s ensitive to the lati t udinal 
distribution of ac tivity (iHowe et all Il999b lAntia et all I2001L 
iHowe et aill2002l) . Therefore the correlation analysis might sug- 
gest that the zonal modes are better correlated with resolved ac- 
tivity index that takes into account the latitudinal dependence of 
solar magnetic activity. To check upon this point, it will be im- 
portant to further extend the correlation analysis to the azimuthal 
components of intermediate - high degree modes investigating 
the degree of correlation with resolved activity index. 



5.3. Conclusions 

Evidences of QBP have been found in the individual low -t 
degree modes as well as in their corresponding azimuthal m- 
components at all levels of solar activity. The strength of the 
QBP signal seems to be modulated by the 1 1 -year envelope, but 
the effect is more significant at low-mid latitudes compared to 
higher ones. This modulation might explain the coupling be- 
tween the strength of the main global dynamo and the appear- 
ance of the QBP mainly over periods coinciding with solar max- 
ima. When the solar cycle is extremely weak, the magnetic flux 
processed by the main dynamo is not strong enough to bring up 
the mid-term periodicity signal in those surface activity indices 
related to the toroidal component of the main dynamo. 

The signatures of QBP have also been detected over the min- 
imum phase of solar cycle 23 with 80% of confidence level. The 
likely persistent nature of the 2-year periodicity in the seismic 
analysis must be investigated, but probably shows a permanent 
and periodic phenomenon that stands in the sub-surface layers 
in addition to a pure presence of a magnetic field. The occur- 
rence of the QBP signal over periods of solar maxima as well as 
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solar minima implies that the origin of the QBP is not coupled 
to the presence of active regions and therefore that the mecha- 
nism is different from the one that brings up to the surface the 
active regions. This finding is in agreement with previous results 
obtain ed by analyzing coronal holes data dVecchio & Carbonel 
l2008h . 

During periods of low activity the QBP signal from low-mid 
latitudes is almost absent, while it appears at higher latitudes. 
This finding might imply that the signal strength increases with 
latitudes and also that the sources of excitation might cover all 
range of latitudes. This feature might be explained within the 
formalism of a second dynamo mecha nism situated near the sub- 
surface layer dBenevolenskaial,[l998a|[bl) . 

The seismic investigation brings new insight on the features 
of solar magnetic activity that might be useful to improve dy- 
namo models. To analyze the signature of the QBP further in- 
vestigation of the sectoral and zonal components of higher de- 
gree modes is strongly desired. We should also carry out a more 
detailed investigation on the latitudinal dependence of the signal 
strength over solar cycle. The findings could be extremely use- 
ful to fully identify the physical mechanism behind it, as at this 
early stage we cannot completely rule out the magnetic Rossby 
instability or other dynamo model s that already take into ac - 
count a further periodic component dFluri & Berdvuginall2004h . 
Moreover we would be able to check whether or not the peri- 
odicity differs in the two h emisphere as already pointed out by 
IVecchio & Carbonel (120081) . 

The end of the solar cycle 23, characterized by its long ex- 
tended minimum and the ongoing solar cycle is a current subject 
of research. These new solar cycle features has re-opened in- 
terests and debates on the dynamo mechanism that rules solar 
magnetic activity. 
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Table 1. The Spearmann correlation coefficient and significance between Activity Indeces and low degree />mode frequency shift 
over different phases of solar activity. 
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Table 2. The same as Table 1 but for the azimuthal component of low degree modes. 
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